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• Remote Sensing of the Water Budget • Coupled Modeling : Diagnostic and_ Prognostic 

Where is the water and how much is there? Where is it going today, where will it go tomorrow? 
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Now, what do we do with it? 
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• Overarching Motivation and The Result 
A Improved Water and Energy Management 
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> Overview collaborative radar rainfall project between Tennessee 
Valley Authority (TVA), VCSI (The Von Braun Center for Science & 
Innovation), N/*|1 a MSFC and ultHuntsville 

• AREPS: ARMOR Rainfall Estimation Processing System 

Demonstration project of real-time radar rainfall using a research radar 

• NREPS: N EXRAD Rainfall Estimation Processing System 

Expansion to Tennessee River Valley using operational WSR-88D 

• Objectives, methodology, some results and validation, operational 
experience and lessons learned. 



NASA Earth Science - Applied Science Program’s ARRA Project: 
Water Supply and Management in California 

• NASA MSFC and UAHuntsville to provide NREPS rainfall input to NASA 
MSFC/USRA distributed hydrological model for soil moisture and 
evapo-transpiration estimation. 


IteMlias, 16r|^E|/futui§ mark 
CA and beyond. 




Objective and Motivation 

> TVA River management and distributed rainfall 
measurements- reducing dependence on rain 
gauges. 

• Provide custom-tailored, radar-based rainfall products 
specific to TVA’s operational river management 
needs (e.g., 6-hr sub-basin mean rainfall) 

• Potential reduction of costs associated with 
maintenance of large TVA rain gauge network 

> Research & Development Operations 








TVA River Forecast Center 

http://www.tva.com/river/flood/center.htnn 
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Knoxville, TN: staffed around the clock and 365 days a year 


River management duties include 

• Issuing forecasts of reservoir levels and water releases at TVA dams 

• Providing hourly generation schedules for TVA hydroelectric projects 

• Providing special notifications to the public during flood events 

• Evaluating cooling water needs for TVA coal-fired and nuclear plants 

• Monitoring water quality conditions below TVA dams (oxygen levels) 

• POC in event of a river system emergency 


Complex, interdependent and sometimes conflicting set of river 
management requirements and responsibilities. 


• Forecasters are busy! 
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Tennessee River Watershed 



emo project using ARMOR 


> TVA River 
Scheduling 
Division 


> 1 12 sub-basins 
^ 1 840 Im 2 

> 189 rain gauges 
maintained by 
TVA 

• Annual costs: 
$6K / gauge 


Tennessee River Basins 


dam and mouth of Tennessee river, 710 Square Miles!^ 





Objective: 

Transition from rain gauge 
(point) estimation paradigm to 
radar (distributed) measurement 


Rain Gauge-Dominated Rain Estimation 


CURRENT River Operations 
Gauge Rainfall Inputs 


Advanced Radar, QPE 
Applications 


FUTURE TVA 


River Ops Rainfall Inputs 


Radar-Dominated ( reduced gauge } Rain Estimation 




Point (gauge) vs. areal (radar) rainfall estimates 


20070706: 24hr Rainfall Accumulation 
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24 Hour Rain Totals July 6, 2007 



Rain Gauge (inches) 


BUT much of the heaviest precipitation missed at individual gauges (this is 
fairly typical)! 

Heterogeneity of rain field presents problems for ptHntrneasurements, even 
for gauge-adjusted radar maps ( ( 



Conventional Radar Estimates of Precipitation 




• Raindrops (often the target) 
are deformed by drag as 
they fall. 

• Become more oblate with 
size (diameter) 


• Radar scans rain/precipitation by sending out a series of microwave pulses along a 1° 
beam while sweeping that beam through a 360° circle. 


• Executed for several elevation angles and repeated every 5-6 minutes. 


• Returned power from rain (called reflectivity- “Z”) is typically measured at only one 
polarization (e.g., horizontal) and related to rainfall via “Z-R” power laws: Z = aR b 




QPE: Problem with conventional radar-rainfall approaches: 
Reflectivity Factor (Z) - Rainfall Rate (R) Relations 


Problem: Numerous rainfall-reflectivity relationships, 
which one is correct? 


Random error up to 100% (instantaneous) can occur. 
Typical space-time smoothing reduces to say 20%-40%. 
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Measurement sensitive to rain drop size distribution, 
presence of hail/ice/snow, and radar calibration. 
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Without care, unacceptable errors/bias can be incurred 
for high resolution hydrological applications. 
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Even gauge corrections are still beholden to gauge 
calibration/error/sample mismatch- a problem at times 
(more often then we would like to admit). 


30 



Sample of current operational relationships: 

Z = 300 R 1 4 - convective rain 
Z = 250 R 1 2 - tropical rain 
Z = 200 R 1 6 - summer stratiform rain 
Z = 130 R 2 0 - winter stratiform (eastern US) 
Z = 75 R 2 0 - winter stratiform (western US) 

How do we correct these issues? 


50 vs 100 mm/hr 
at 50 dBZ over a 
valid range of 
observed DSDs!! 



Z (dBZ) 




Enter Dual-Polarimetric Radar 


• Approach: Provide measurements 
in both horizontal and vertical 
polarizations (dual-polarimetric). 

•Advantages: Better description of 
various particle types/shapes in a 
given volume 

•Determine size distribution- more 
accurate rain rates (improved QPE) 

•Rain drop shape related to size 

•Hydrometeor ID (hail vs. rain) and 
non-meteorological scatterers 
(clutter!) 

•Consistent calibration 




We need the 
measurement in H 
and V directions! 







Enter Polarimetric Radar 


And more variables 

1 . Reflectivity factor Z at horizontal (Z h ) or vertical (Z v ) 
polarization [Conventional radar measure] 

- Measure of drop size and concentration; 

• Most sensitive to drop SIZE (D 6 ) 

2. Differential reflectivity Z DR (Z h /Z v ) 

- Measure of median drop diameter— > SIZE/SHAPE 

- Useful for rain / hail / snow discrimination— ► SIZE/SHAPE 

3. Specific propagation differential phase K DP E(k h - k v ) 

- Measure of water content and drop size— > NUMBER/SHAPE 

- Immune to radar calibration, attenuation, partial beam blockage 

4. Correlation coefficient p hv 

- Indicator of mixed precipitation — > SHAPE/PHASE/CANTING 
(Depolarization) 

- Useful for identifying non-meteorological scatterers too! 



TRANSMIT 

Horizontal Pulse 
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Advantages: Better description of various particle types/shapes in a given volume 

• Determine size distribution- more accurate rain rates (improved QPE) 

• Hydrometeor ID and non-meteorological scatterers (clutter!) 

•Consistent calibration 
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Dual-Polarimetric Radar: Improved QPE through 
improved description of particles 


Particle Size-Controlled 


Smalle r ZDR Larger ZDR 


Small Drops Large Drops 


Rain is Oblate 


Small Drops (1 mm) 

Large Drops (> 4 mm) 

Axis ratio decreases 
with increasing size- 
more oblate 


Insects 


Hail/Graupel Rain 


Radar “sees 
an effective 
sphere 


Hail/Graupel 


Smaller KDP Larger KDP 


Number Concentration Controlled 


Axis ratio ~ 1 

Tumbling and lower dielectric strength makes 
hail look like spheres 

Unless they start to melt. . . 

Melting Hail/Graupel 
(Toroid or ice core; 

looks like a huge drop) ^jg ratj() < 1 


Smaller# 


Larger # 


Larger ZDR Smaller ZDR 


Large Drops Small Drops 


Smaller KDP Larger KDP 







Advanced Radar for 
Meteorological & Operational Research 



> Location: 

• Huntsville International Airport, Huntsville, AL 
(Altitude 206m) 

> C-band dual-polarimetric Doppler radar 

> Simultaneous transmit and receive of H, V 

> Variables: Z, V, W, ZDR, 0 DP , p hv 

> Operations: 

• 24-hrs a day / 7 days 

• Rain volumetric scans at least every 5-min 
(tilts: 0.7 ,1.5 ,2.0 ) 

• Also used in research mode (e.g., RHIs, full 
volumes, vertically pointing scans) 

> Routine calibration: 

• Receiver calibrations 



ARMOR Rainfall Estimation 
Processing System (AREPS) 
























Distance North fkm'l from ARMO 


AREPS Coverage 

> 1 00 km from 
ARMOR (HSV) 

>11 sub-basins 

> 42 rain gauges 




and State* 



AREPS Distributed Rainfall Products 




'TV A Detain: 6-I3F Rain stats [inches! ? IOC Jena P.odiTis from AFJTGP. 
Dare/ Time |EfTC|: 2UU31Z 1U_Q55957 

Esn# mean Haw H in S-Dev * kP^JL 

3603 0.434 0.375 0,139 0,529 100,000 

3602 U.769 1.523 0.29® 2.050 100 . UUU 

3701 0.351 0.616 0.006 0.S5O 99.112 

3605 0 r 7?6 1.969 0.13P 1,046 100,000 

3702 0.672 U.904 0.009 1.167 100 . 000 

3601 1.265 £.034 0.46B 3,096 103,000 

3604 1.246 2.4DB 0.196 0,522 100-000 

3703 0.9S7 1.727 0.332 1.342 100 . 000 

3504 0.239 1.13B 0.051 0,320 103,000 

3501 0.120 0.349 0.000 0-434 37.500 

2 SOL 1.040 1.466 0.609 4.445 100 . 000 


http://vortex.nsstc.uah.edu/armor/webimaqe 


> Rainfall products created every 5-min: 

• 1 -hr and 6-hr basin/sub-basin rainfall 
statistics (mean, max, min, etc) 

• Rainfall at critical locations (e.g., dams) 

• rainfall accumulation images (1 -hr, 6-hr) 

> Text files transmitted every hour to TVA 

• Contain previous hour’s rainfall 

• used as input by inflow model input 
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Alabama 


Distance East (km) from ARMOR 


With River Basins and States 


6-hour accumulation statistics 


3703 


Distance East (km) from ARMOR 


With River Basins and Rain Gauges 




Verification: Point Comparisons 
AREPS vs. TVA rain gauges 
(October 2007 - June 2008) 


Before Calibration Correction 


After Correction 
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Bias = -17% (-1 .80 mm) 
Error = 18% 


75 


T 


Bias = -10% (-0.99 mm) 
Error = 12% 


Radar Rainfall Estimate Improved 
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Verification: Sub-basins 
AREPS vs. rain gauge-derived areal mean 
(JanuaHf 2008 - July 2008) 
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Gauge derived accum. (mm) 


gaugemetwork with respect to 
stfbbasih boundaries) 

Demonstration successful. . ..but what about the rest of the watershed? 


Expansion of radar rainfall maps 
across entire TN River watershed 



Numbar of 
radars 

3 



NEXRAD coverage across TVA watershed grid 
(Number of radars within 200 km of a gridpoint) 
(Yellow and red gauges are owned by TVA; red = critical) 
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With River Basins and Rain Gauges 


Distance East (km) x rc- 

^ ' -> Sufficient radar coverage 





NEXRAD Rainfall Estimation 


Processing System (NREPS) 



4) Merge radars 


2) QC/rain rate algorithm 


NCEP 


3) Gridding 2 km 2 
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WSR-88D Level II 
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NREPS Processing 


> VCP-based sweep selection (e.g., 1 st SUR and Doppler scans of 
each elevation)) | 

> Beam blockage correction for occultation 

> Clutter mitigation (N EXRAD ORPG 3) 

• Notch width filter 

• Sun strobe removal 


> Non-precipitation mitigation (Steiner and Smith 2002) 

• Vertical extent and structure 

• Spatial variability 


> Additional checks using VCP mode and RUC analysis melting level 

• Clear-air VCP + 0 C height > 100 m -> no precipitation 

• Clear-air VCP + 0 C height < 100 m -> snow 


> Rainfall from individual radar interpolated to 2 km x 2 km Cartesian 
grids using NCAR reorder at 1 .0, 1 .5 and 2.0 km ARL for TVA (more 

se lowest available “qood” rainfall 
lndividu< 



Beam Blockage Correction 


> DEM vs ray path 
model -> regions of 
occultation 

> Visibility correction 

factor (Germann et al. 2006) 

f = (1 -blockage)' 1 

> Adjust reflectivity field 
for blockage of 5% to 

90% (Lang et al. 2009) 


Terrain within 200 km of Radar (m) 
1-1:25 (-1i24 -1|23 — =m — =321 — =020 — =1 



Beam-blockage fraction 
Mi25 (-1S4 -123 — =WS — =0j21 — =120 — =0=1 



o.oofop 

KDAX 0.5 elevation 

Worst case scenario: Azimuth: S3" 
BBF at 200 km: 0.51853533 
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Non-precipitation mitigation 
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Steiner and Smith (2002) 
Reduce 


Decision tree methodology 

performed in polar coordinates 

1 ECHOtop -> echoes separated 
from precifflHHMHHHHHH 

2 SPINchange -> echoes 
embedded in precip 

3. vertGRAD -> limits echo removal 
around edges of storm cells 


Computationally fast and 
efficient 

ifficulties with widespread 
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NREPS Rain Rate Algorithm 



R = R (Z h ) / SLWE 

SLWE = snow /liquid 
( SLWE = 8 for TN Valley) 













NREPS Performance 



> Error = 20% 

> Bias = -7% 

> R 2 = 0.93 

> NREPS has been 
implemented by TVA 

• used routinely in river 
scheduling operations 


> TVA rain gauge removal 

. Gauges selected after 




Coverage: Radar vs Rain gauge 



• Intense isolated rainfall in 
west-central TN (sub-basin 
3911) 

•Amounts: (6pm-midnight) 

- Gauge = 38.1 mm (1.5 in) 

- Radar = 125 mm (4.9 in) 


• Which was right ? 

- radar beam < melting level 

- max reflectivity < 50 dBZ 

- stationary thunderstorm 

- no rain gauges in 3911 


^ r . Advantage : Radar 






“1000-year flood” in Middle and West Tennessee 
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Downtown Nashville TN flooding, 
especially near Cumberland River, 
which crested at 52 ft there, and 
tributaries. 


Gauges with Theissen 
polygons to-sub-basin 6-hr 
mean to basin mean 

NREPS grids to 
sub-basin 6-hr means to 
basin mean. 


6-hour periods (Daily quarters) 


12.0 


10.0 


id.u 


T 2T 


0.05 


20100503 0000 UTC N EXRAD 24— HR RAINFALL ACCUMULATION 


-300 -200 -100 0 100 200 300 

Distance East (km) from NEXRAD 


1 00 


-100 


-200 


200 




With Sub- basins 


2 May 2010: NREPS 24-hr Rainfall 
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Source: Mr. Adam Cissna. TVA 


6 Hour Valley Totals 

Tennessee River Basin: April 30-3 May 2010 
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In many sub-basins, there is good agreement between TVA-gauge and NREPS 
6-HR cumulative rainfall. 


Sub-basin 6-HR Rainfall 


Sub-basin Cumulative Event Rainfall 


Pickwick Savannah 



Pickwick Savannah Cumulative 
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Perryville Johnsonville 



Perryville Johnsonville Cumulative 
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TVA Cumulative ^^—Radar Cumulative 


Figures are courtesy ofIVIr. Adam Cissna, TVA 






More good agreement... 


Sub-basin 6-HR Rainfall 


Savannah Perryville 
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Sub-basin Cumulative Event Rainfall 
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lures are courtesy ofMr. Adam Cissna, TVA 






And some more (including one with 


Sub-basin 6-HR Rainfall 


Cheatham to Clarksville 



•TVA 


•Radar 


Scottsboro Guntersville 
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Sub-basin Cumulative Event Rainfall 


Cheatham Cumulative 
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Scottsboro Guntersville Cumulative 
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ires are courtesy ofIVIr. Adam Cissna, TVA 








But, as seen in the Basin mean, there was a tendency for the radar < gauge 
rainfall in some sub-basins, especially in “2 nd wave” of heavy rain on May 2 nd . 


Sub-basin 6-HR Rainfall Sub-basin Cumulative Event Rainfall 



Cbcrrtesy Mr. Adam Cissna, TVA 







Some more sub-basins with radar < 


Sub-basin 6-HR Rainfall 


Duck above Manchester 



■TVA Radar 


Normandy 



■TVA Radar 


in 2 nd peak... 


Sub-basin Cumulative Event Rainfall 
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Duck above Manchester 



Normandy Cumulative 



lures are courtesy ofIVIr. Adam Cissna, TVA 





And there are almost always a couple of exceptions to the rule! 
Radar > gauge, especially in the 1 st wave of heavy rain on May 1 st . 


Sub-basin 6-HR Rainfall 
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Sub-basin Cumulative Event Rainfall 



Courtesy Mr. Adam Cissna, TVA 
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TVA-VCSI Project Practical Experience 


> 


> 


> 


Detailed developer (user) knowledge of user requirements, priorities 
and computing environment (tool capabilities and limitations) is key. 

• With diverse backgrounds, takes time and committed project team to develop 
productive communication and working relationship to solve complex issues. 

• Active user participation in system testing and improvement is very helpful. 

In an operational hydrological setting like the TVA River Forecast 
Center, rainfall data availability is the primary priority. 



• WSR-88D data pathways (Internet and LDM stability, back-up), 

• Software design (exception handling, stability, quality assurance) 

• Hardware design (memory, CPU and storage specifications to worst case 
scenarios; redundant systems) 

Data accuracy is important but secondary. 

• Because of river management complexity and model 
uncertainty, it is difficult for TVA to provide 6-hr sub- 
basin rainfall accuracy guidelines. 



Peak Discharge 
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stream flow gauge Sata in real tirneS Model output is 
tuned to measured river response. 
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California Water Resources Project 



> Funded by NASA through 
American Reinvestment 
& Recovery Act 

• Through June 201 1 

> Develop integrated tool to 
assist water management 
within the San Joaquin 
River Valley 

. Radar precip estimates 

. Distributed hydro model 

. Snowfall measurements 

• Surface temp / moisture 



Challenges 


> Residual non-precip 
contamination 

> Nature of precipitation 

• Shallow, non-bright brand 

• Snow vs rain 

> Terrain and remaining 
occultation 

• Non-standard refraction 

> Beam height vs melting 
level 


Map of lowest usable beam height in Central CA 





1-HR Rainfall Amount (inches) 


Radar vs rain gauge comparison 
during mixed-precip events 


Jan 29-30 2010 

Radar vs Gauge (GID: 383 Petersburg, TN) 
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Future NREPS plans 


> Vertical profile of reflectivity (bright band mitigation, ice 
vs. rain, rain rate vs. beam height) 

> Explore more sophisticated precipitation type (C/S) 
identification methodology 

• Other categories for Coastal California - shallow, non- 
convective, non-bright band ? 

• z-R’s for the appropriate regime and precipitation type 

> For CA ARRA NREPS: Integration with NASA distributed 
hydrologic model for soil moisture and evapo- 
transpiration estimates 


> Rank tn tha fiitnra — As W.RR-RRD natwnrk is nnlarizarl 



